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Abstract 

Besides the loop contribution and the coupling changing, heavy Higgs decaying to 125GeV Higgs 
is another important effect in Higgs property measurement in the models that have extend Higgs 
sector. Such heavy Higgs decay effect contain more than one sources, such as A ^ Z/i, H /i/i, 
A 2 Aih^ and Hs H 1 H 2 . Some of them will be promising in NMSSM Higgs sector searching. 
Higgs data fitting is performed with and without this heavy Higgs decay effect being included, 
and the final result demonstrate such effect is very important in Higgs data fitting if the signal 
efficiency is large. 
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I. INTRODUCTION 


With a 125GeV Higgs-like scalar found by the LHc[l]j2], the last cornerstone of the 
Standard Model(SM) have been confirmed. The measurement of the scalar’s couplings to 

nn 

vector boson and fermions are compatible with the SM predictionsEven so, due to 
some difficulty like the gauge hierarchy, dark matter, or vacuum stability, it is commonly 
believed that SM is just an effective theory in low energy scale and New Physic(NP) will 
emerge in some higher energy scale. While the LHC haven’t find any unambiguons NP 
signals, so a main work at present is using the low energy observables, especially the Higgs 
data obtained these years, to constrain the NP models’ parameter space. A convenient and 
popular method to perform the Higgs data constraints is through the Higgs signal strength 
fitting. Higgs signal strength is the ratio of observed Higgs signal number to SM predicted 
Higgs signal number, and Higgs signal number is the product of Higgs production crosssec¬ 
tion, branching ratio, and signal efficiency. If the NP model have a light particle that can 
change the Higgs effective coupling to photon and gluon in the loop level, or have a scalar 
that can mix with the 125GeV SM-like Higgs, or lead to a non-SM Higgs decay channel, then 
the Higgs signal strength will change, then one can use the observed Higgs signal strength 
to exclude the NP points that leading to a too large or too small signal strength. In snch 
sitnation, the topology of the Higgs process are just the same as the SM Higgs process. It 
means, you don’t need to worry about the signal efficiency, because signal efficiency have 
been cancelled completely in the ratio (the signal efficiency in the nnmerator and denomina¬ 
tor are equalled). Such argument is true when the non-SM Higgs production crosssection 
is negligible compared to the SM Higgs production crosssection. But, if the non-SM Higgs 
prodnction crosssection and SM Higgs production crosssection are on the same level, then 
the effective coupling and branching ratio are not enough for the Higgs signal strength fit¬ 
ting. The signal mixing of the SM produced Higgs and non-SM produced Higgs should be 
included. While such a possibility have been analysed in Minimal Supersymmetry Stan¬ 
dard Model(MSSM) (Higgs and neutralino associated production) and two Higgs doublet 
mode 1(2HDM)3(SM -like Higgs are decayed from heavier scalars). In this paper we will 
analyse the non-SM Higgs prodnction in the Next to Minimal Snpersymmetry Standard 
Model(NMSSM) and its influence on Higgs data fitting. Now let me explain why we choose 
NMSSM. 
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Among various New Physic models, Supersymmetry (SUSY) is a very attractive one due 
to its elegant solution to hierarchy problem and other theoretical and phenomenological 
merit like providing dark matter candidate, unifying the gauge coupling, or stabilizing the 
vacuum in high energy scale. The most studied and popular low energy scale SUSY model 
is MSSM. In MSSM, there are 5 Higgs states(2 CP-even neutral scalars, 1 CP-odd neutral 
scalar and a pair of charged scalar) after spontaneous symmetry breaking, and the lighter 
CP-even scalar can be regarded as the 125GeV scalar found in LHC. But the up limit of 


the lighter CP-even scalar mass in MSSM is about ISOGeVlTl, so the 125GeV Higgs will 

n 

leads to some degree of fine-tuning[8||. The /r term in MSSM super potential also cause some 
theoretical puzzle. While in NMSSM, fine-tuning can be alleviated by an additional singlet 


which can enhance the tree level Higgs massj^. At the same time, the puzzling /i term can 
be explained by a dynamical origin. More paper related to naturalness, see . On the other 
hand, if the non-SM Higgs production process is heavier scalar’s decay, then there always 
need a small tan /jQjfl]. It is partly because small tan [5 will lead to a small branching ratio 
to bottom quark pair. Besides, if the heavy scalar is heavier than 2 times of top quark mass, 
then the branching ratio to top pair will be dominant and the branching ratio to 125GeV 
Higgs will be strongly supressed. Such an small tan (3 and light non-SM Higgs scenario is 
unfavoured in MSSM but can be easily obtained in NMSSM(l2l| |l8l| Ibecause the Higgs mass 
enhancement in tree-level and the mixing with the singlet). So the NMSSM Higgs sector is 
very hopeful in non-SM Higgs production. 

In this paper, we will analyse the production and decay of the NMSSM Higgs sector in 
detail, show how large this non-SM Higgs production crosssection can be, and estimate the 
impact of this non-SM Higgs production on Higgs data fitting. In next section, we will briefly 
review the NMSSM Higgs sector. In section 3, we will detailedly describe our scan method. 
Section 4 is used to show the results, and we analyse the Higgs data fitting including non-SM 
Higgs production in section 5. Section 6 is the conclusion. 


II. NMSSM HIGGS SECTOR 


The NMSSM Higgs sector contain two Higgs doublets and a Higgs singlet jl^: 


H„ = 




H,= 


H7 


,A. 


( 1 ) 
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Without regard to domain wall problem, the Z3-invariant NMSSM superpotential and cor¬ 


responding Higgs soft breaking terms are 




1 


W = u‘^KQB^ - - BKLHa + XSH^Ha + 


( 2 ) 

( 3 ) 


-C^oft = mj,jHu\^ + mj,jHd\^ + mi,\S\^ + [XAxSHuHd+-tiA^S^ + h.c.]. 

After spontaneous symmetry breaking(SSB), Higgs helds will get vacuum expected values: 



Vd 


v2 \ 0 


. {S) = 


( 4 ) 


with V = \/v‘^ + v'^ = 24:6GeV. 

Using three minimum conditions, and m| can be represented by other param¬ 

eters, then the independent parameters used to describe NMSSM Higgs sector are: A, k, 
/i = Aus, tan/3 = u„/ud, A\^ A^. Rotating the unphysical goldstone state away, the mass 
matrix of pseudoscalar is: 

An, 


Ml,, 


= Ml, (Ml = + KV.)) 


Mpi2 = - 3AkuJ —, 

Z Vg 


Mp22 


= - {M\ sin 2/3 + SXkvA — SkVsA,^/\ p2. 
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( 5 ) 

( 6 ) 
( 7 ) 


Our work is concentrate on the Higgs sector, so it will be more convenient to replace A\ and 
A^ by Ma = Mpii and Mp = Mp 22 in NMSSM Higgs sector description. Ma and Mp are 
the mass scale of heavier doublet and singlet separately. But unlike MSSM, Ma and Mp 
are just two approximate mass scale in NMSSM, not physical observables. 

Then the scalar mass matrix is: 

1 


Min = Mi + {Ml - -X^v^)siM2/3, 


M‘si2 


= -i(A/|-iAV)sm4/3, 


, ^9 1 / , ,9 . 2 . 21 1 9x ^ COS 2 B 

Msi3 = -^(^ASm2/3 +AkuJ-, 

Z Vg 

M‘s22 — Af| cos^ 2/3-|--(An)^ sin^ 2/3, 

2 


M; 


523 = ^(2A^u^ - M\siT? 2/3 - Xpv^ sin2/3) —, 

Z Vg 


1 ^2 2 
Ms 33 = -^M\sm^ 2(3—-\-2i^v‘l-\-KVgA^/s/^ —-^Xkv^ sm2(3. 


( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 
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The charged Higgs mass is: 


= Ml + MS, - i(At.f. (14) 

Where Mz and Mw are the mass of Z boson and W boson separately. The scalar mass 
eigenstates are denoted by H 2 , and H 3 with increasing mass. The pseudoscalar mass 
eigenstates are denoted by Ai and A 2 with A 2 is heavier than Ai. Charged Higgs is denoted 
by CH. 

A analytical study can be performed in MSSM Higgs sector (tree-level) if Ma is heavy or 
the tan/3 is large than 5, because the Higgs spectrum and the Higgs mixing angle can be 
expressed by or While in NMSSM, Ma and Mp maybe not too larger or smaller 

than V, and tan (5 tend to be small, so such an expansion is impossible in most situation. In 
addition, NMSSM Higgs sector is described by 6 parameters, and the physical observables’ 
dependence on these 6 parameters are complicated. So we will scan the parameter space that 
leading to what interests us, and use a lot of result to show the property of this parameter 
space. In next section we will describe the scan method in detail. 


III. SCAN METHOD 


We use NMSSMTools4.5. 1 jf7| t o scan the parameter space. Package NMHDECAY|18| and 
NMSDECAy[i^, based on HDECAY ^ and SDECAyI^I, are used to calculate the decay widths 
and branching ratios of Higgs and sparticle. The NMSSM particle mass spectrum (include 
Higgs mass), mixing angle, and conplings are also calculated by NMSSMTools4.5.1. 


A. Scan range 


The value range of the 6 parameters mentioned before can be confine d by a series of 
conditions snch like vacuum stability or perturbativity in high energy scale (l6|]. Combined 
with the argument in the introduction, we set our parameters range as below: 

• In order to avoid Landan Pole before the Great Unification scale, A and k, can’t be too 
large: 

0.0<A,k<0.7. (15) 
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• As we mentioned in the introduction, small tan/? is favoured by large Higgs decay 
branching ratio, so: 

1.0 < tan/3 < 5.0. (16) 


• In MSSM, the lower limit of Ma is about 300GeV(except in the alignment region|22l|). 
While in NMSSM, due to the mixing with SM singlet, the lower limit of Ma can reach 
200GeV. Besides, we don’t want too heavy non-SM Higgs, so: 


200GeV <Ma< 500GeV. 


(17) 


• While the singlet mass Mp can vary in a larger range: 


OGeV < Mp< 500Gel/. 


( 18 ) 


• We will focus on Higgs sector, so the stop sector is mainly used to tune the Higgs 
mass. We don’t want a light stop to muss our analysis: 

700GeV < Mq 3 = Mu 3 < 2TeV, -3TeV < A* < 3TeV. (19) 

• /i correspond to Higgsino mass. One may let fi become large to avoid non-SM Higgs 
decaying to Higgsino. But we find a large fi always leading to a too heavy Ma that 
almost touch the 2mt threshold (this can be understood by the pseudoscalar mass 
matrix(5)). So fi in our scenario is not too large : 

lOOGeV <ii< bOOGeV. (20) 

• All the other soft parameters will be set to 2TeV, and all the A-terms(except At, Ax, 
A^) are equal to zero. 


B. Constraints 


Both the HI and H2 could be the 125GeV SM-like Higgs. So in this work, our discussion 
is performed in two scenarios: 

• HI scenario: Hi is the 125GeV SM-like Higgs. 

• H2 scenario: H 2 is the 125GeV SM-like Higgs. 
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NMSSMTools have included lots of low-energy observables from B-physics, LEP, Tevatron, 
and LHC to constrain the parameter space. Besides those traditional low-energy observables 
constraints, we describe some related constraints used in our work. 


There is only one Higgs in 122GeV-128GeV(we have excluded the degenerate Higgs 
scenario in this work(22t]), and its Higgs signal strength should be consistent with Higgs 
data in 2cr range. 


HiggsBounds-4.1.2 [2^ is used to constrain the Higgs sector. 


Non-SM Higgs searching through H/A —> rr in CMS 25] 26] and ATLAS 28], and 
through H —>• AA in CMS 27] have been included in NMSSMTools. In addition, the 
observed limit given in the non-SM Higgs search through A ^ Zh and H ^ hh is 
used to constrain the parameter space also. 


The thermal relic density of the 
upper bound measured by Planck 


ightest neutralino is required to be lower than the 
3nl ] , and the rescaled proton scattering cross section 


need to satisfy the direct detection bound from LUXjMl- The calculation is performed 


by MicroOmega 


32j. 


The lightest sfermion should be heavier than SOOGeV. 


Our scan is a multi-step scan. First, we scan the entire parameter space quickly, then 
we intensively scan the part that will leading to significant result. Finally we sum all the 
partial results together. In next section we will show our scan results and analyse them in 
detail. 


IV. SCAN RESULTS 

In this section we will show our scan results and discuss why these plots are like that. 
We will show the production crosssection of heavy scalars first, then analyse their branching 
ratio, finally the combination of crosssection and branching ratio will tell us how many 
125GeV SM-like Higgs are produced from heavy scalars’ decay. 
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FIG. 1: Left plot is the total production crosssection of H 2 and H 3 at 8TeV LHC, right plot is the 
total production crosssection of Ai, A 2 , and CH at 8TeV LHC. These results are obtained in HI 
scenario(iLi is the 125GeV SM-like Higgs). 


A. The Production Crosssection of Scalars 


The gluon fusion process gg ^ H and gg ^ A are calculated by HIGLU(33t] in NNLO 


While for vector boson fusion , top pair associated production, and associated production 
with vector boson, we just use the results inj^ . And due to the tan (3 in our scenario is small, 
contribution from bottom quark fusion can be safely neglected. For the production cross 
section of charged Higgs, we use the results inj^]. Combined with the reduced couplings 
calculated by NMSSMTools, we can get the production crosssection of all the 6 Higgs states. 

We want to explore how many 125GeV Higgs can decay from heavy scalars, so we will 
show the production crosssection of the scalars that are heavier than 125GeV only. FigHJis 
the total production crosssection of H 2 , H 3 , Ai, A 2 , and CH at 8TeV LHC in HI scenario, 
and FigElis the total production crosssection of H 3 , Ai, A 2 , and CH at 8TeV LHC in H2 


scenario. Compared 


20pb at 8TeV LHC 


o the total production crosssection of the 125GeV SM Higgs(about 
3J]), the doublet dominant H3 and A2 still have a signihcant production 


crosssection. The larger production crosssection of pseudoscalar is due to the fermion loop 
factor in Higgs-gluon-gluon effective coupling, and the peak around 2mt can be explained 
by the mixing effect of scalar and top quarkoniumj^. So, if the H^ and A 2 branching ratio 
decaying to 125GeV Higgs are large, then such heavy Higgs decay effect may be important 
in Higgs signal analysis. The singlet dominant H 2 and Ai may also contribute to this 

















FIG. 2: Left plot is the total production crosssection of H^, at 8TeV LHC, right plot is the to¬ 
tal production crosssection of Ai, A 2 , and CH sd, 8TeV LHC. These results are obtained in H2 
scenario(iL 2 is the 125GeV SM-like Higgs). 


effect if they mixing with the doublet adequately. CH can contribute to this effect through 
CH hLF(we will use h to denote the 125GeV SM-like Higgs in below). But the production 
crosssection of CH is too small(because CH can’t be produced through gluon-gluon fusion, 
its production will be suppressed by the bottom quark PDF and the associated produced 
top quark in the hnal state), and CH will mainly decay to top bottom(depend on tan /5)(35|. 
So we don’t expect a large contribution from CH. 


B. The Branching Ratios of Scalars 

Higgs SSB is the origin of SM fermion and vector boson mass. So bottom pair, top pair, 
W boson pair and Z boson will be the main decay product of the neutral Higgs. But the 
pseudoscalar can’t couple to gauge boson in lowest order, so the pseudoscalar branching 
ratio to WW and ZZ can be ignored. The Higgsino and Singlino in our scenario are not too 
heavy, so the branching ratio to and can not be ignored. The sum of the branching 
ratio to and will be denoted as BRS: 


3 


BRS{H/A) = + 

(21) 

3 

BRS(CH) = Y.BR(x°xt) 

(22) 


1=1 
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In order to illustrate the heavy Higgs decay effect more obviously, we deffne a quantity 
”Numh”. Numh is the average number of 125GeV SM-like Higgs produced from the heavy 
scalar decay. For example, if Numh{A 2 ) is 0.7, then 10000 A 2 will produce 7000 125GeV 
SM-like Higgs on average. The cascade decay process like A 2 —?• H 2 Z —?• hhZ also need to 
be included in Numh calculation. A detailed describe of Numh calculation can be found 


in jo]!^]. We emphasize here that Numh only include decay topology containing at least 
one 125GeV SM-like Higgs, but sometimes the heavy Higgs will tend to decay to non-SM 
Higgs. In such situation Numh will be small even the other branching ratios are small. 

FigEland ffgjUare the braching ratios of heavy Higgs in HI and H2 scenario respectively. 
A main property illustrated by these plots is that signiffcant Numh can only appear below 
the 2mt threshold, as we mentioned before. Light singlino and higgsino occupy a large part 
of decay product, especially in the H2 scenario. H 3 always have a possibility to get a large 
Numh before its mass reaching 2mt threshold. Due to a large production crosssection, A 2 
can contribute largely to the heavy Higgs decay effect in HI scenario. While A 2 in H2 
scenario tend to decay to non-SM Higgs, so the contribution from A 2 will be sub-dominant 
in H2 scenario. CH’s branching ratio is mainly occupied by BR(tb) and BRS. So combined 
with the production crosssection, contribution from CH will be very small. 

We have production and decay of the NMSSM Higgs sector have been shown clearly, next 
we will show how much 125GeV Higgs can be produced from the heavy scalar’s decay. 
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FIG. 3: The main branching ratios of heavy Higgs in HI scenario. 
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FIG. 4: The main branching ratios of heavy Higgs in H2 scenario. 


C. Heavy Higgs Decay Effect in 125GeV Higgs production 


Result can be simply obtained by multiplying production crosssection and Numh. But 
before giving the final result, we will show the specific channel that have a significant con¬ 
tribution. Besides the decay channel, the production process are separated also. We find 
all the important contribution are from H 3 and A 2 in both scenarios. FigiS] and FigEl are 
the results. Those processes could be search channel of NMSSM Higgs sector at LHC run 


H 


38l | ^ . But we need to notice that such a significant crosssection is sensitive to the mixing 


between non-SM and SM Higgs, and the non-SM Higgs mass scale. The mixing and mass 


scale allowed 
measurement 




current Higgs data will be excluded largely by future precise Higgs property 


40(] ( 1 ^ ^ . So before we perform the NMSSM Higgs sector searching in future 
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FIG. 5: Channels that contribute significantly in HI scenario. The gluon-fusion crosssection is 
obtained ar 8TeV LHC. 



FIG. 6: Channels that contribute significantly in H2 scenario. The gluon-fusion crosssection is 
obtained ar 8TeV LHC. 

collider, these channels might be not promising anymore under updated constrain. 

FigJT] is the sum of all production and decay channel for all the scalars that are heavier 
than 125GeV Higgs. As we mentioned in introduction, this heavy Higgs decay effect is 
strongly suppressed by large tan /3. This effect is much significant in HI scenario that in H2 
scenario. In next section, we will analyse how this effect change Higgs data fitting. 
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FIG. 7: The sum of all production and decay channel of all the scalars that are heavier than 
125GeV Higgs. Left plot is obtained in HI scenario and right plot is obtained in H2 scenario. 

V. THE IMPACT ON HIGGS DATA FITTING 


In the final section, we will show how the heavy Higgs decay effect affect the Higgs data 
fitting. Considering the uncertainty in the gluon-gluon fusion process calculation is larger 
than 15% and the production crosssection of SM Higgs in 8TeV LHC is about 20pbj^, we 
will focus on HI scenario only. Let’s briefly review the Higgs data fitting first. 

The data used to fit is the Higgs signal strength p which is the ratio of observed signal 
event number to the SM expected number. While the NP models may lead to a p that not 
equal to 1, so the observed /r can be used to constrain the NP models, /i can be calculated 
by production crosssection, branching ratio, and signal efficiency: 




a{X)BR{h Y)e 

o'SM{X)BRsM{h — 7 - Y)esM 


(23) 


Where X is production process, Y is the final state, e and esM are the signal efficiency in NP 
model and SM respectively. If the non-SM Higgs production process can be ignored, then e 
and esM can be ignored too, because they will be just the same|43||. 

In our HI scenario, the heavy Higgs decay effect can not be ignored (larger than 5pb), 
so the calculation of p will be different. In such situation, signal efficiency can not be 
ignored again, we need to know how the non-SM Higgs production process feed each of 
the signal region measured by the experiment. Dividing the numerator and denominator 
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simultaneously, [i will be: 


l^{X,Y) = 


a{X)BR{h ^ y) + a{HD)BR{h Y) 


e 

^SM 


(24) 


<7sM{X)BRsM{h —>• Y) 

'''HU'' denote heavy Higgs decay, and we will use HD for short in the remaining part. It’s 
obviously that the core problem in here is to decide the ratio of efficiency so a detailed 
collider simulation is needed. But we do not prepare to perform a collider simulation in 
this work. We will just suppose several different and show how the Higgs data fitting 
change. The decision of —^ by detaifed simufation wilf be our future work. By the way. 


the report refeased recently |4j] show a discrepancy between the experimental observation 
and theoretical prediction in Higgs process, so a precise simulation seems not easy. On the 


other hand, our fitting process will not be performed to the exclusive signal region like 


451 . 


We will use the combined result of ZZ and 77 channel given by the Fig.2 of jsj to perform 
the fitting. Then can be calculated simply by: 

2 _ - ^ 77 )^ I - i^ZzY 

X AO 


A 2 

77 




(25) 

Where the hatted number and A are the central value and error bar given by Fig .2 of j^. 

Next we will show the distribution in tan/3 — Ma plane in three scenarios: HD is 
ignored; —^ is 0.4; —^ is 0.8. But, could the ratio of efficiency —^ be as farge as 0.8 or 
0.4? There are two reason for it: 


• The mass spfitting between Ma 2 and Mh + Ma^, or Mh^ and Mh + Mhj^/Mh 2 are not 
too large, so the differential distribution of the 125GeV Higgs produced from HD will 
not be too different from the 125GeV Higgs produced from SM processes. 

• The main decay product of Ai, Hi, and H 2 are bottom pair or neutralino, but the 
description of the incfusive signal of ZZ and 77 channel (section 5.1 of [? ] and section 
V of [? ]) do not mention requirement on b-jet or missing energy. 

Due to these two reasons, we think HD will feed the Higgs searching signal effectively. (We 
need to emphasize again that our work just give a rational assumption, the precise value of 
—^ can only be obtained by collider simulation.) 

FigJHl is the comparison of these three scenarios. In order the illustrate the HD effect 
more cfearly, we pick the points with HD effect larger than l.Opb. The color pattern show 
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FIG. 8: The tan /3 — Ma plane 

a clear changing in Higgs signal strength fitting. Dne to BR{h —)■ ZZ) and BR{h —)■ 77 ) 
suppression (because hbb will be enhanced by mixing with heavy doublet), and the mixing 
with singlet, the signal strength of h in NMSSM always tend to be smaller than the SM 
expected value. So the HD's feeding always lead to a smaller or say higher confidence 
level. 

VI. CONCLUSION 

Higgs signal strength fitting always ignore the signal efficiency. But this method is right 
only if the non-SM Higgs production is very small. In this work, we find A —>■ Zh, H —)■ hh, 
A 2 — 7 - Alh, and H 1 H 2 can contribute largely to the 125GeV SM-like Higgs production. 

The corresponding total crosssection can be large than fi.Opb. Then we discnss how such 
a non-SM Higgs prodnction process will change the Higgs data fitting. By a reasonable 
assnmption to the signal efficiency, the fitting result change largely by the feeding of this 
non-SM Higgs production process. The precise estimation of this effect need detailed collider 
simulation which beyond the scope of this work. 
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